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Abstract. In this paper we consider the ignition and structure of a reacting compressible mixing layer 
using finite rate chemistry lying between two streams of reactants with different freestream speeds and tem- 
peratures. Numerical integration of the governing equations show that the structure of the reacting flow can 
be quite complicated depending on the magnitude of the Zeldovich number. In particular, for sufficiently 
large Zeldovich number, the three regimes first described by Linan and Crespo (1976); ie., ignition, 
deflagration, and diffusion flame, occur in supersonic as well as in subsonic flows. An analysis of both the 
ignition and diffusion flame regimes is presented using a combination of large Zeldovich number asymptot- 
ics and numerics. This allows us to analyze the behavior of these regimes as a function of the parameters of 
the problem. For the ignition regime, a well defined ignition point will always exist provided the adiabatic 
flame temperature is greater than either freestream temperature. One important result is that at supersonic 
speeds ignition occurs far downstream from the plate and, as the flow is accelerated to hypersonic speeds, 
ignition is exponentially delayed. For the diffusion flame regime, the location of the flame changes 
significantly with changes in the equivalence ratio and the Schmidt numbers. 
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1. Introduction. Quite recently it has been realized that an understanding of the structure and the sta- 
bility characteristics of compressible, particularly supersonic, reacting mixing layers is extremely important 
in the context of the scramjet engine. We have begun a systematic study (Jackson and Grosch, 1990) of the 
inviscid spatial stability problem of these flows with the combustion zone modeled by the flame sheet, as 
described by Jackson and Hussaini (1988). We found that the addition of combustion in the form of a flame 
sheet had important, and complex, effects on the flow stability. In addition, we also showed that with 
sufficient heat in the flame sheet the flow can undergo a transition from convective to absolute instability, 
even without the presence of a region of reversed flow in the layer. There are of course situations in which 
the flame sheet model may not be a sufficiently accurate representation of the combustion process. This 
could be important in the calculation of the stability characteristics of the reacting flow. Therefore in this 
paper we return to the problem of determining the structure of the laminar diffusion flame in a compressi- 
ble, particularly supersonic, mixing layer in using finite rate chemistry. 

In a seminal paper, Linan and Crespo (1976) examined the structure of a diffusion flame in the 
unsteady mixing of two half spaces of fuel and oxidizer. By using a combination of large activation energy 
asymptotics and numerics they showed that three laminar regimes exist within the flow, the ignition, 
deflagration, and diffusion regimes. These same regimes were also found by Jackson and Hussaini (1988) in 
their study of combustion occuring between two parallel supersonic streams with nearly the same speed and 
temperature. These three regimes are shown in Figure 1. The ignition regime is a region where the com- 
bustible gases mix until, at some finite distance downstream of the plate, a thermal explosion occurs and the 
gas is ignited. The second regime is the deflagration region. After ignition, a pair of well-defined 
deflagration waves (or "premixed flamelets") emerge according to classical thermal explosion theory. These 
waves penetrate the mixing layer until all of the deficient reactant is consumed. Just downstream of the 
deflagration wave, a diffusion flame regime exists where the mixing process is governed by diffusion in the 
direction normal to the flame. This regime covers the bulk of the overall flow field. Collectively, the igni- 
tion and deflagration regimes can occupy a very small portion of the flow field, especially at supersonic 
speeds. These regimes are highly nonlinear and hence are not easily amenable to analysis. If the incoming 
flow were turbulent rather than laminar, we surmise there would be a collection of initially discrete 
premixed flamelets. Full numerical solutions, based on the compressible Navier Stokes equations, appear to 
miss the ignition and the deflagration regimes, and hence, the premixed flamelet is not captured. Perhaps 
careful fine-tuning of the numerics will uncover these regimes. 

The spatial development of the reacting compressible mixing layer is considered here without making 
the assumption of infinite Damkohler number. If the Damkohler number were to be taken to be infinite, the 
flow and combustion fields could be described by similarity solutions. With a finite Damkohler number, the 
flow and combustion fields evolve in the downstream direction and only approach the similarity solutions of 
the flame sheet at infinite distance from the plate. The mean fields with finite Damkohler number will, in a 
future paper, be used in stability calculations. 

In the next section we formulate the problem. Numerical solutions to the full equations are then 
given in section 3. Sections 4 and 5 present an analysis of the ignition and diffusion flame regimes using a 
combination of large activation energy (and hence Zeldovich number) asymptotics and numerics. Conclu- 
sions are then given in the final section. 
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2. Formulation of the Mean Flow Equations. The nondimensional equations governing the steady 
two dimensional flow of a compressible, reacting mixing layer with zero pressure gradient lying between 
streams of reactants with different speeds and temperatures (Figure 1) are given by (Williams, 1985) 


(p U) x + <pV) y =0, 1 = p7\ (la,b) 

p {UU x +VU y ) = (iiU y ) y , (lc) 

p (UT X + VT y ) = Pr~ l diT y ) y + (y - \)M 2 \iU y 2 + pn, (Id) 

p(UF^ x + VF hy ) = Scf i ( i iF hy ) y -^a, j = 1,2, (le) 

n = D pF l F 2 e' ZtlT . (10 


In these equations, the x axis is along the direction of flow; the y axis is normal to the flow; U and V are 
the velocity components in the x and y directions, respectively; p is the density; T is the temperature; and 
Fi and F 2 are the mass fractions, with the reaction assumed to be irreversible and of the Arrhenius type. 
The viscosity p is assumed to be a function of temperature. The nondimensional parameters appearing 
above are the Prandtl number Pr , the Schmidt number Scj = Pr Lej for species j where Lej is the Lewis 
number for species j, the Mach number M = U^l , the Zeldovich number Ze = EIR with E the 
dimensional activation energy and R the universal gas constant, D the Damkohler number defined as the 
ratio of the characteristic diffusion time scale to the characteristic chemical time scale, P the heat release 
per unit mass fraction of the reactant, (3, a parameter involving stoichiometry, and finally y the specific- 
heats ratio. The equations were nondimensional i zed by the freestream values T w p^, £/«, F } for the tem- 
perature, density, velocities and mass fractions, respectively. Lengths are referred to some characteristic 
length scale of the flow. 

The boundary conditions consistent with (1) are 

T = U = Fj = 1, F 2 = 0 at x = 0, y > 0, and x > 0, y °o, (2a) 

p' 

T = p r , U = py < 1, F, = 0, F 2 = = (0 2 / 3i) 4» _1 of x = 0, y < 0, and x > 0, y -> - °=, (2b) 

r 

where <J> is the equivalence ratio defined as the ratio of the mass fraction F x in the fast stream to the mass 
fraction F 2 in the slow stream divided by the ratio of their molecular weights times their stoichiometric 
coefficients. If <j> = 1, the mixture is said to be stoichiometric, if <)> > 1 it is Fj rich, and if <|> < 1, it is F x 
lean. If p r is less than one, the slow gas is relatively cold compared to the fast stream, and if p r is greater 
than one it is relatively hot. For what follows we take the stoichiometric parameter to be unity. 

The mean flow equations (1) are first transformed into the incompressible form by means of the 
Howarth-Dorodnitzyn transformation 
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y = 


y 



y 

V = pV + f/j p, dy. 


(3) 


Next we transform to 


T| -r\ 0 


y_ 

~ 2<x 


(4) 


which is the similarity variable for the chemically frozen heat conduction problem, and T | 0 corresponds to a 
shift in the origin. For the case of both streams being supersonic, T| 0 is determined uniquely by a compati- 
bility condition found by matching the pressure across the mixing layer, resulting in 


1 

Vm 2 - l 


V{x,°o) = 


-1 P(/ 

Pf 1 - 1 Pr 


V (*,—), 


which reduces to that given by Ting (1959) for (3 r = 1. We note here that if the freestream velocities were 
subsonic, then the compatibility condition would be trivially satisfied, and thus r\ 0 would remain indeter- 
minate (Klemp and Acrivos, 1972). Under these transformations, with U =/'(t|) and V = (j\f -f)Hx, 
equation (1) become 


/** + 2 //" = 0, (5a) 

Axf'T x - Pr~ x T" -2fT'-(y- \) = Ax^D F x F 1 e- Zerr , (5b) 

4 xf'Fj' X - Scf ] Fj - 2 fFj = -Ax D F i F 2 e~ Ze>T , (5c) 

where the primes indicates partial differentiation with respect to the similarity variable T|, and we have 

assumed the linear viscosity law \i = T (for Chapman’s linear law \x — CT y the constant C can be scaled 

out by rescaling rj and / appropriately, but must be borne in mind when transforming the variables back to 
their dimensional forms). In terms of the transformed variables, the boundary conditions are 

T = /' = F x = 1, F 2 = 0 at x = 0, ti > 0, and x > 0, i) -> «>, (6a) 

T = p r ,/'= P(/, F { = 0, F 2 = 4> _1 at x = 0, r\ < 0, and x > 0, ti -> -°°- (6b) 

Finally, we note that the Damkohler number D can be scaled out of the equations by a rescaling of the x 
coordinate. Thus, the Damkohler number is not a true parameter of the above system. However, rescaling x 
by the Damkohler number is not particularly useful. Since D is typically exponentially large the rescaled 
coordinate will be exponentially stretched, and this is not desirable from a numerical viewpoint. 

3. Numerical Solutions. The system (5) with boundary conditions (6) must be solved numerically. 
Noting that / is only a function of rj, we see that (5b, c) are parabolic in x and thus can be solved by a 
marching procedure in x. We do this by a Crank-Nicolson procedure. This leads to an implicit nonlinear 
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system of equations in T] which must be solved at each x step. These are solved by an iterative scheme 
using the boundary conditions in q as given in (6). 

In Figure 2a we plot the maximum temperature in T| versus x for Ze - 10, M = $ u = 0, 
Pr = P = <t> = 1* The adiabatic flame temperature is 1.5 for these conditions (see equation (50)), and is larger 
than either freestream temperature. There is a smooth and gradual transition from the inert solution at 
x = 0 to a diffusion flame. In Figure 2b we show similar results but with Ze = 30 and all other parameters 
the same. Note the sudden transition between the inert solution and the diffusion flame solution occurring 
in a narrow range in x around x = 3. This initial region defines the ignition regime, and hence the flame 
location. Similar results are shown in Figure 3 but for M - 2 with an adiabatic flame temperature of 1.7. 
From this figure it is apparent that the ignition regime also exists in supersonic flows. As the Mach number 
increases there is a corresponding increase in the inert temperature at x = 0 due to viscous heating, thus 
lessening the relative effect of combustion on the overall temperature field. 

In Figure 4 we plot the temperature and mass fraction profiles as a function of T| at various x loca- 
tions for the same conditions as in Figure 2b. In Figure 4a the rapid rise in the temperature over a narrow 
range of x should be noted; this indicates ignition. Note also that the temperature peak shifts and that the 
profile is asymmetric due to the asymmetry in the Lock velocity profile. Figure 4b shows the profile of the 
F i mass fraction. At the smaller values of x t F\ decreases gradually from 1 at t| = + oo to 0 at q = -°o, 
with F t varying most rapidly over the range -3 < T| <1. This shows that there is a diffusion of F \ from 
the i\ > 0 region into the T] < 0 region with F\ = 0 only at q = — ®o. As * is increased the rale of change of 
F x with T] over this zone increases. For larger x one can see that F j drops to zero for r| just below zero 
showing the presence of a diffusion flame in which all of F\ is consumed. But note that there is a small 
secondary maximum in the F\ distribution in — 3 < T) < — 1 . This shows the presence of a "premixed 
flamelet" in this premixed region. The distribution of F 2 in r\ at various x is shown in Figure 4c and is the 
converse of that of F x . The only exception is in the secondary maxima, showing that the presence of the 
f premixed flamelet" is much more pronounced. Again the asymmetry of the secondary maxima is due to the 
asymmetry of the velocity profile. 

The existence of the "premixed flamelcts" and the diffusion flame is shown quite clearly in Figure 5a, 
where we plot the loci in the (x,T|) plane of the maxima of Q, the chemical production term, also for the 
data of Figure 2b. As shown in this figure, the position of the maximum decreases from about T| = -0.2 at 
x = 0 to nearly T] = -0.6 at x ~ 2.9. The maximum value of Q increases along this curve. At x = 2.9 
ignition occurs and two maxima appear giving rise to the "premixed flamelets". Beyond the ignition point 
the premixed flamelets move outwards in the shear layer until all of the deficient reactant is consumed. The 
appearance of the third maxima just behind the ignition point marks the appearance of the diffusion flame. 
As x is further increased, the diffusion flame becomes dominant and, as x — > the diffusion flame thins 
and approaches a flame sheet characterized by local chemical equilibrium (Jackson & Grosch, 1990). Note 
that the diffusion flame location is constant in q as x increases. 

In Figure 5b we show similar results for the same values of the parameters except that p 7 = 0.5 and 
P = 1.5. This case corresponds to unequal freestream temperatures. The adiabatic flame temperature is 1.5 
(again see (50)), and is larger than either freestream temperature. Note that the ignition point has moved 
into the region of higher freestream temperature, and the location of the diffusion flame is unchanged while 
that of the premixed flames has changed. In Figure 5c we show results for p 7 = 0.5 and P = 0.4. The 
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adiabatic flame temperature in this case is 0.95, and is smaller than the freestream temperature at + °°. In 
contrast to the two previous figures, there is no well defined ignition point; the premixed flame merges 
smoothly into the diffusion flame whose location is unchanged. In addition, die premixed flamelets are 
absent. Finally, we note that the behaviour described in Figure 5 also occurs for Mach numbers greater than 
zero. 


In the following two sections, an analysis of the ignition regime and a brief discussion of the diffusion 
flame regime is given using a combination of large activation energy (and hence Zeldovich number) asymp- 
totics and numerics. 

4, Ignition Regime. To investigate the above behaviour analytically for large Zeldovich numbers, we 
take Lej = 1, Pr = 1, resulting in Scj = 1. This allows us to consider linear combinations of (5b) and (5c) 
which eliminates the source terms. The solutions of the resulting equations, using the boundary conditions 
(6), are given by 


T + PF, = |} r + (1 - p r + + :L y-M 2 { 1 - M 2 ^! (7) 

t + pf 2 = p r + P0- 1 + (i - p r - pr 1 )'? + 1 y Lw2 ( 1 - (8) 


where 


vj/ = 


(T-PvVO-M 

(1 + erf T | )/2 


0 < p y < 1 
Ptf = 1- 


( 9 ) 


Thus F i and F 2 can be eliminated from the energy equation (5b), yielding 

4xf’T x -T"-2f r'-(y- l^a-py) 2 ^') 2 

= [Pr + (1 - Pr + P)^ + Zy-MW ~ Pr/) 2 ^0 - 'P) - T] 

x [pr + Pr 1 + (1 - Pr - P<t>"’)'P + - Py^'PO - *¥) - 7']<r z " T (10) 

for the temperature, subject to the boundary conditions 

T = 1 at x = 0, r| > 0, and x > 0, q -> °°» (1 la) 

T ~ at x = 0, rj < 0, and x > 0, q — > -oo. (lib) 

In the absence of chemistry, equation (10) reduces to 


4xf’T x -T" -If f-(y- l)M 2 (\ - P^O?') 2 = 0, 


( 12 ) 
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which possesses the inert similarity solution 

T, = T, (n) = 1 - (1 - Pr)(l -'¥) + 1 - pt,) 2 ^ - 40, (13) 

which is also a solution of (10) for small x. As x increases, more of the combustible mixes until, at some 
finite distance downstream of the plate, a thermal explosion occurs characterized by significant departure 
from the inert. 

To analyze the ignition process, we expand about the inert by defining T } 

T^ZeiT-Tj), (14) 

and take the asymptotic limit Ze — » «>, Substituting into (10), we see that the left hand side becomes 
0(Ze ! ) while the right hand side is 0(D e *), where T* = max(Tj ). In order to balance these terms, 
D must therefore have the form 


D=-^-e 2 *' T - 

P5Ze 


(15) 


where 5 will be chosen in the course of the analysis. Substituting (14) and (15) into (10) yields the follow- 
ing equation for Tj 


4* T X x - T\ - 2/ T\ =4* g|j[P 'F-Ze-’r,] 

x [p<r'(l - V P) - Ze“ ! Fi ] exp 


T\ + 2e(r y - T.) 
T. (Tj + Ze~ l T{) 


(16) 


which must be solved subject to the boundary conditions 

T x = 0 at x=0 and at T| — > ±<x>. (17) 

Note that the right hand side of (16) is exponentially small except where |7) - T* | = 0(Ze~ l ). There are 
now several cases to consider, depending upon the magnitude of the parameters {J r , and M . Each case 
is described below. 

4.1. |l-|5 r |«l, |l-p e/ |«l > and M=0( 1). 

This case corresponds to the simple but important case of ignition for nearly equal freestream temperatures 
and velocities, and the analysis is presented by Jackson and Hussaini (1988). For this reason, we only briefly 
review the results here. We begin by setting 


$T=Ze( 1 - P T ), % = (1 - Py), 


(18) 
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with and M fixed and 0(1). In particular, we see that 


/ = r\ + 0{Ze~ v2 ), V = 'V o +0 (Ze~ m ), '¥ t> = | (1 + erfr\). 


( 19 ) 


The inert solution (13) is expanded as 


T, = 1 + Ze" 1 [-P T (1 - yj + ^y^A/ 2 ^ T, (1 - '¥„)] + 0(Ze~ 2 ), 


( 20 ) 


which yields T+ = 1. Substituting into (16), together with the choice 5=1, yields the following equation 
for the disturbance temperature 


4xT x -2x\T[ -Tf = 4x V F 0 (1 - x ¥ 0 ) exp 


T x - p r (1 - VJ + *y± A/ 2 ^ W 0 (1 - V 0 ) 


( 21 ) 


subject to the boundary conditions (17). This equation was solved numerically by a Crank-Nicolson implicit 
scheme. Asx increases, the solution for T { becomes unbounded at some finite location (;c*,T]*) downstream 
of the plate, indicating thermal runaway. The results are presented in Figure 6, where we plot x* and r[* as 
a function of M for various values of P 7 . Note that for a fixed p^, increasing M causes the ignition 
point to move closer to the origin. For a fixed M increasing (or decreasing) from zero (i.e., cooling 
(or heating) the gas from below) causes r\* to move into the region of higher freestream temperature. In 
addition, x * increases as the gas is cooled below and decreases as the gas is heated below. 

As noted above, as p r increases or decreases from zero, the ignition point moves into the region of 
higher freestream temperature. This case corresponds to different freestream temperatures, and can be 
analyzed in the limit 1 I — > °°. We shall assume p r < 1 (p r > 0) and hence the ignition point is located in 
the region t| >> 1. The analysis for the case p r > 1 (p 7 < 0) is similar. Following the analysis of Linan and 
Crespo (1976), we set 


j5 r (l-^) = z, 0<z<°°, (22) 

and using the asymptotic expansion of *¥ 0 , we see that 

H 2 = In / 2 Vi) - In T| - In z. (23) 

The reaction zone is located where z =0(1), and to a first approximation we get 

r|. 2 = In (jij- / 2 Vi) - In T}» , (24) 

where T|. defines the ignition location. Since the chemistry is confined to a thin region about T). , we shift to 
the reaction zone coordinate 


1 

2t|. 


In z. 


T1 =T1 


(25) 
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Substitution into (16), together with 7* = 1, and neglecting terms of 0(qr 2 ), yields a canonical equation for 
T\. Repeating the analysis for [3/ > 1 yields a canonical equation for T x of exactly the same form, the 
differences being in the definition of the parameters, the flame now being located in the slow stream where 
the gas is initially hotter, and T, = p 7 . This canonical equation, originally given by Linan and Crespo 
(1976) in a different context, is 


-x { 1 z 2 V„ = (z - Ay)e' 1 ' - * 


(26) 


where 


*i =' 


x / (Pr ri* 2 ) Pr < 1 

x / (~ P 7 T| 2 ) Pr > 1. 


A = 


<K1-Pr)/P 

(Pr-D/P 


Pr < 1 
Pr > 1. 


(27) 


This equation must be solved subject to the boundary conditions 


V(0) = 0, v.(~) = 0. (28) 

The first boundary condition is a restatement of (17), white the second boundary condition comes from 
matching with the chemically frozen region q = 0(1). Note that x is now a parameter of the equation. 

Linan and Crespo have shown that for A < 1, is doubled valued for x l <x u , single valued at 
xi = x\ m> and does not exist for x { > x\ 9 . Near x\ % we note that 5\(//3xi becomes large showing that the x 
partial derivative can no longer be neglected in (16). Therefore X \ 9 may be considered as the ignition point. 
For A > 1, \|/ is single valued for all Jt] indicating that there is no well defined ignition point. In this case 
the premixed flame merges smoothly with the diffusion flame and there are no premixed flamelets. This 
behaviour is shown numerically in Figure 5c. In Figure 7 we plot the ignition point X\ m versus A. As A 

approaches one from below, the ignition point moves off to infinity, while for A > 1 there is no well defined 
ignition point. 

4.2. |l-Prl<<l. Pt/ < 1, and M« 1. 

This case corresponds to nearly equal freestream temperatures in a shear flow at small Mach numbers, and 
we therefore set 


= Ze (1 - p 7 ), M=JZ£M, (29) 

with P r and M fixed and 0(1). The inert solution (13) is now expanded as 

T, = l + Ze" 1 [-pr(l- v f0+ 1 ^ 1 W 2 (l-p (/ ) 2 ‘R(l-‘P)] + O(Ze- 2 ), (30) 

again yielding T* - 1. Here, 'F is given by (9a). Substituting into (16), again with the choice 5=1, yields 
the equation for the disturbance temperature 
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4x f T\ - If T[ - T" = 4x 4>(1 - »F) exp 


r, - p T (1 - 40 + ^A/’d - Pt/^'Pd - 40 


, (31) 


subject to the boundary conditions (17). This equation was again solved by a Crank-Nicolson scheme. As x 
increases, the solution for T x becomes unbounded at some finite location (x*,q*) downstream of the plate, 
indicating thermal runaway. Some results are presented in Figure 8, where we plot x * and q* as a function 
of P(/ for various values of M and p T = 0. We first note that for = 1, the ignition point is independent 
of the Mach number. For a fixed M , decreasing P^ from one causes the ignition point to move back 
toward the plate and into the slower moving gas. At fixed p y < 1, an increase in M causes the ignition 
point to shift toward the origin. 

Increasing or decreasing P r has the same effect on the ignition location as in the previous case. To 
see this, we first note that the asymptotic expansion of the Lock profile (9a) for large q is given by 


/' = 


1 - A (q + B ) 1 exp (- (q + B ) 2 ) 

Pi/ +A m + BT'expi-Pum + B) 2 ) 


q -> +°o 
q -> -oo, 


(32) 


where A and B are known constants. Following the analysis of the previous section, we again obtain the 
canonical equation (26) for the temperature perturbation, with x x now given by 




x / (p r (q* + B ) 2 ) Pr < 1 

*/(-Mq. +£) 2 ) P t > h 


(33) 


and A as before. 


4 . 3 . p r >0, 0<Pt,<l, and M > 0. 

For Mach numbers greater than zero, p^ < 1 and any value of p r , the inert solution (13) has a maximum at 
T{ = 0; i.e.. 




1-Pr + X y i M 2 d~Pf/) 2 d-24') 


= 0 . 


(34) 


If 'P =0, then the location of the reaction zone will occur in the outer region |q | >> 1 where the gas is 
initially hotter. Since the analysis for this case will be similar to that of the previous two cases, we will not 
present it here. If the second term of (34) is zero, then the reaction zone will occur in the mixing layer, and 
is given by the implicit relation 


^ *'P(Tl.)= y + 


1 ~ Pr 

(7- l)Af 2 (l -Pt,) 2 ' 


(35) 
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The corresponding maximum value of 7} is given by 


T* = max (7/) = 1 - 


\ ~ Pr y~ 1 


Af 2 (l-M 2 + 


(1 - Pr) 2 

2(y- 1) A/ 2 (1 - Pu) 2 


(36) 


Substituting (34-36) into (16) shows that the right hand side of (16) is exponentially small except in a nar- 
row region where |7) - T» | ~ 0 (Ze~ ] ). Outside this region the equation is essentially frozen, and is given 
by 


4-* T l Jt = Tj” + 2/ T[ , (37) 

valid on either side of the reaction zone. To study the structure of the reaction zone we define the new 
coordinate 


r\ = T). + 8z, 


(38) 


where 5 is the thickness of the reaction zone and is not yet specified. The temperature and velocity are also 
expanded as 


T i = T u (x) + 87' u (x,ti) + 


*P = »P. + 8z + 


(39) 


Substitution into (16) yields 


d}T 

ii = -4 X *P. (1 - 'POexp 

dz 2 V 




2,xv'\2,2 


TS 


(40) 


where we have chosen 8 = Ze . Note that (40) can be integrated once to yield 


x) J -J=-exp[-b 2 z 2 ]dz + c,(x). 


(41) 


where 


d( x) = 4xW. (1 - 'P.)exp[T u /T. 2 ]^, (42a) 

0 

b 2 = M 2 ( 1 - p t/ ) 2 ep;) 2 7’.- 2 (42b) 

Evaluating (41) at z = ±°o yields the matching conditions for (37). Specifically, we now solve (37) subject 
to the boundary and jump conditions 



- 11 - 


T 1 (x,°°) = 0, 7’ 1 (j:,ti. + ) = T u (x), 


dr , 

dTl 


(X,T1. + )=C1(JC)-C/(X), 


(43) 


Ti(x,-°°) = 0, i ( j: , T( . ) = ri.(x). 


dT, 

dil 


(jc. 11. ) = Ci(x), 


with the initial conditions 


(44) 


r^o.Ti)^ o, q*q*. (45) 

Equation (37), valid on either side of the reaction zone, can be solved numerically by marching forward in 
x. The unknown functions T u (x) and c x (x) are chosen at each x location so as to satisfy the boundary con- 
ditions at T| = ±oo. Thermal runaway will occur at (x*,q*)» thus defining the flame location. In Figure 9 we 
show T Xm and Cj as a function of x for p 7 = 0.5, Pf/ = 0, and M = 2. Note that as x approaches x *, T Xm 
becomes unbounded indicating thermal runaway. In Figure 10 we plot the ignition location x * and q. as a 
function of the Mach number for p y = 0 and various values of p 7 . As the Mach number becomes large x* 
moves off to +°o while q* tends to zero for any p r . This shows that in the supersonic region ignition 
occurs far downstream of the plate and, as the flow is accelerated to hypersonic speeds, ignition is exponen- 
tially delayed. For unequal freestream temperatures (P r * 1), as the Mach number goes to zero then r \ * 
moves out of the mixing layer into the region of higher freestream temperature. When this occurs this 
analysis is no longer appropriate, and one must consider the case of *¥ = 0 in (34). 

5. Diffusion Flame Regime. Typically, a diffusion flame is characterized by a chemical reaction time 
that is much smaller than a characteristic diffusion time. Chemical reactions then occur in a narrow zone 
between the fuel and the oxidizer, where the concentrations of both reactants are very small. Mathemati- 
cally, the assumption of very fast chemical reaction rates leads to the limit of infinite Damkohler number 
which reduces the diffusion flame to a flame sheet (i.e., local chemical equilibrium). This assumption 
significantly reduces the complexity of the problem since it eliminates the analysis associated with the 
chemical kinetics. For many flows, the assumption of local chemical equilibrium adequately predicts the 
location and the shape of the diffusion flame (Buckmaster and Ludford, 1982; Williams, 1985). Since the 
flame sheet model has a simple similarity solution and hence is amenable to analysis, we give its structure 
below. For finite values of the Damkohler number, equations (5b, c) must be solved directly using the 
numerical scheme mentioned above. 

Assuming unit Lewis and Prandll numbers and a linear viscosity law, the profile associated with the 
flame sheet model is given by (Jackson and Hussaini, 1988, Jackson and Grosch, 1990) 

= i — (i + <t> -1 )(i - '?), r 2 = 0, (46) 

T = P r + P* -1 + (1 - Pr - P^ 1 )^ + 1 y lA/2 ( 1 " Pi/) 2 ^ 1 - '*0. ( 47 > 


valid for q > i \ f , and 
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Fi = 0, F 2 = qr ] -0 +<tr 1 )'P, (48) 

T = p r + (1 - p r + p)¥ + A/ 2 (l - P t/ ) 2 'P(1 - 40, (49) 

valid for i] < 7]y- . Here, gives the location of the flame sheet where both the reactants vanish, and T 
takes the adiabatic flame value 


T f = Pr + (1 - Pr + P)^/ + ^ M 2 (l ~ P u) 2y V f (1 - '*'/)• (50) 

The flame location is given by the implicit relation 

^/*'p ( % ) = T77 (5i) 

which is independent of P^ and M . This is only valid for the case of a linear viscosity law where the 
momentum equation decouples from the energy equation. For a more general viscosity law, the flame loca- 
tion will depend on p r and M as well. When using (3) to transform back to the physical coordinates, Lhe 
diffusion flame location in physical space will depend on p^, <(>, p r and M , independently of the viscosity- 
temperature relation. The internal structure of the diffusion flame, which appears as a discontinuity on the T| 
scale, is described in in Jackson and Hussaini (1988). 

The above assumption of unit Lewis and Prandtl numbers led to a closed-form analytical solution. 
Since the Prandtl number is about 0.72 and the Lewis numbers are not one, the Schmidt numbers will in 
general not be one. To determine the influence of the Schmidt numbers on the location of the diffusion 


flame, we begin by considering (5b,c) in the limit of infinite Damkohler number 

Pr~ l T" + 2/T +(y- 1) M 2 (f") 2 = 0, (52) 

Scf' Fj + 2 fFj = 0, (53) 

valid on either side of the flame sheet. The appropriate boundary and jump conditions are given by 

T - 1, ^1 = 1, F 2 = 0 as ti (54a) 

[ Pr~' T + P Scf' F-] = 0, [T] = 0, Fi=F 2 = 0, at T\=r\ f , (54b) 

T = Pt-, Fj = 0, F 2 =<() _1 as i) — ^ (54c) 

where [ ] denoles the jump ( )^ + - ( ) at r\j , Integrating the equations (53) for the mass fractions yields 

Fi = 1 -cJcP') Sc, rfTi. F 2 = 0, (55) 
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valid for r| > r|/ , and 


F, = 0, F 2 = <r‘ -c 2 j &f 2 dn, 
valid for tj < t\j . Using the boundary condition F ] = F 2 = 0 at T)y yields 

i r 1 

c ' = ^ * c *-~ f 

J C¥'f'dr\ j Q¥) SCl dr\ 


( 56 ) 


(57) 


The location of the flame sheet 71/ , where both reactants vanish simultaneously, is given by the implicit 
relation 


* 1 / 

Q>Sc 2 j m SCl dl\ = Sc l C¥}f 2 ' Sc 


J QV'f l d% 

’V 


(58) 


obtained from the jump condition (54b) on the derivatives. In Figure 1 1 we plot the flame location as a 
function of Sc j for various cases. In Figure 11a we show the variation with Sc 2 . As Sc 2 increases the flame 
location moves into the slower moving stream. In Figure lib we show the variation with the equivalence 
ratio. As 4> increases the flame location also moves into the slower moving stream. Finally, in Figure 11c we 
show the variation with (3^. As P v increases the flame moves into the faster stream. 

6. Conclusions. In this paper we have considered the ignition and structure of a reacting compressible 
mixing layer using finite rate chemistry lying between two streams of reactants with different freestream 
speeds and temperatures. Numerical integration of the governing equations show that the structure of the 
reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. In particular, 
for sufficiently large Zeldovich number, the three regimes first described by Linan and Crespo (1976); ie., 
ignition, deflagration, and diffusion flame, occur in supersonic as well as in subsonic flows. Our numerics 
pick up the premixed flamelets as clearly shown by the distribution of the production term. Again, depend- 
ing on the magnitude of the Zeldovich number, there can be a gradual or sudden transition from an inert 
solution to a diffusion flame. 

An analysis of both the ignition and diffusion flame regimes was presented using a combination of 
large Zeldovich number asymptotics and numerics. This allowed us to analyze the behavior of these regimes 
as a function of the parameters of the problem. For the ignition regime, a well defined ignition point will 
always exist provided the adiabatic flame temperature is greater than either freestream temperature. One 
important result is that at supersonic speeds ignition occurs far downstream from the plate and, as the flow 
is accelerated to hypersonic speeds, ignition is exponentially delayed. For the diffusion flame regime, the 
location of the flame changes significantly with changes in the equivalence ratio and the Schmidt numbers. 

There are substantial changes in the temperature and mass fraction distribution as the mixing layer 
evolves downstream. These are expected to have a considerable effect on the stability characteristics of the 
flow. The changes, and particularly the rate of change, in the mean profiles were shown to be very sensitive 
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to the value of the Zeldovich number. Therefore we expect that the stability characteristics will also be 
quite sensitive to the value of the Zeldovich number. The stability calculations are underway and will be 
reported at a later date. 
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Figure 1. Schematic showing the reacting mixing la 
greater than T x and T 2 . 
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Figure 5. Plot of the loci in the (jc , tj) plane of the maxima of Q for M = Py = 0, 
4> = 1, Ze = 30, and (a) p 7 = p = 1; (b) P r = 0.5, P = 1.5; and (c) p r = 0.5, P = 0.4. 
Here, D = tye 7 * /(P Ze). 
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